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Many  problems  in  acoustics  contain  a  small  dimensionless  parameter  e,  and  it  is  useful,  both 
conceptually  and  from  the  point  of  view  of  numerical  computation,  to  seek  a  solution  in  the  form  of  a 
perturbation  series.  In  the  simplest  cases  the  series  would  proceed  by  integral  powers  of  e, 
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where,  in  a  typical  application,  might  be  an  acoustical  potential  or  pressure  and  e  a  frequency  parameter, 
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perturbation  problem,  that  is,  one  in  which  no  single  series  like  the  one  quoted  wilt  vt  valid  both  in  the 
near  field  (where  boundary  data  are  specified  and  where  surface  loading  may  be  of  interest)  and  in  the  far 
field  (where  the  signal  directivity  and  level  are  required).  Separate  aeries  must  be  developed  describing  the 
near  and  far  fields,  but  neither  can  be  completely  constructed  independently  of  the  other  because  each 
series  lacks  sufficient  boundary  data  for  its  unique  determination. 

This  report  describes  how  the  method  of  Matched  Asymptotic  Expansions  (MAE)  can  be  used  safely 
and  systematically  (1)  to  indicate  the  appropriate  form  taken  by  the  inner  (near  field)  and  outer  (wave 
field)  series  and  (2)  to  determine  all  unknown  functions  and  constants  appealing  in  both  series  by 
“matching”  the  series  according  to  a  clear-cut  rule.  These  points  are  illustrated  by  detailed  study  of 
several  very  simple  problems  in  low-frequency  acoustic  scattering  problems  which  serve  to  demonstrate 
that  physical  arguments  are  unreliable  in  these  problems  and  that  they  are  no  substitute  for  the 
unambiguous  matching  rule.  Two-dimensional  scattering  problems  are  used  to  introduce  logarithmic 
gauge  functions;  it  is  shown  that  the  matching  rule  can  easily  accommodate  these  functions  and  moreover, 
that  insistence  upon  satisfaction  of  the  matching  rule  can  in  some  cases  be  used  to  greatly  improve  the 
rapidity  of  convergence  of  series  involving  logarithmic  functions.  The  report  emphasizes  the  very 
widespread  applicability  of  the  MAE  method  to  problems  in  classical  and  modem,  linear  and  nonlinear 
acoustics  and  related  fields. 
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ABSTRACT 


Many  problems  in  acoustics  contain  a  small  dimensionless  parameter 
e,  and  it  is  useful,  both  conceptually  and  from  the  point  of  view  of 
numerical  computation,  to  seek  a  solution  in  the  form  of  a  perturbation 
series.  In  the  simplest  cases  the  series  would  proceed  by  integral  powers 
of  e, 

0  =  <t>0  +e0t  +e202  +  •••■ 

where,  in  a  typical  application,  <p  might  be  an  acoustic  potential  or 
pressure  and  e  a  frequency  parameter,  or  Helmholtz  number.  In  most 
cases,  however,  the  problem  of  determining  the  functions  0j  is  a  singular 
perturbation  problem,  that  is,  one  in  which  no  single  series  like  the  one 
quoted  will  be  valid  both  in  the  near  field  (where  boundary  data  are 
specified  and  where  surface  loading  may  be  of  interest)  and  in  the  far 
field  (where  the  signal  directivity  and  level  are  required).  Separate  series 
must  be  developed  describing  the  near  and  far  fields,  but  neither  can  be 
completely  constructed  independently  of  the  other  because  each  series 
lacks  sufficient  boundary  data  for  its  unique  determination. 

'>  This  report  describes  how  the  method  of  Matched  Asymptotic 
Expansions  (MAE)  can  be  used  safely  and  systematically  (1)  to  indicate 
the  appropriate  form  taken  by  the  inner  (near  field)  and  outer  (wave 
field)  series  and  (2)  to  determine  all  unknown  functions  and  constants 
appearing  in  both  series  by  “matching”  the  series  according  to  a  clear- 
cut  rule.  These  points  are  illustrated  by  detailed  study  of  several  very 
simple  problems  in  low-frequency  acoustic  scattering  problems  which 
serve  to  demonstrate  that  physical  arguments  are  unreliable  in  these 
problems  and  that  they  are  no  substitute  for  the  unambiguous  matching 
rule.  Two-dimensional  scattering  problems  are  used  to  introduce 
logarithmic  gauge  functions;  it  is  shown  that  the  matching  rule  can  easily 
accommodate  these  functions  and  moreover,  that  insistence  upon  satis¬ 
faction  of  the  matching  rule  can  in  some  cases  be  used  to  greatly  improve 
the  rapidity  of  convergence  of  series  involving  logarithmic  functions.  The 
report  emphasizes  the  very  widespread  applicability  of  the  MAE  method 
to  problems  in  classical  and  modem,  linear  and  nonlinear  acoustics  and 
related  fields.  c 
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1.  INTRODUCTION 


Singular  perturbation  techniques  have,  over  the  past  25  years,  been  extensively  applied 
to  problems  in  fluid  dynamics,  the  field  in  which  many  of  these  techniques  were  first  intro¬ 
duced  and  in  which  most  of  their  development  has  taken  place.  The  literature  relating  to  the 
method  of  Matched  Asymptotic  Expansions  (MAE)  alone  is  vast,  running  to  six  books  [  1  -6| , 
several  dozen  review  articles  and  more  than  three  thousand  research  papers,  as  at  mid  1976. 
Despite  the  prolixity  of  these  papers  in  many  branches  of  steady  and  unsteady  fluid  mechanics, 
it  is  indeed  surprising  that  modern  perturbation  methods  should  have  made  so  little  impact  on 
one  particular  branch  of  unsteady  fluid  mechanics  -  Acoustics.  It  is  even  more  surprising  when 
one  recalls  the  work  of  the  great  1 9th  Century  workers  in  acoustics  -  Helmholtz,  Kelvin, 

Stokes,  Rayleigh  -  much  of  which  has  many  ideas  in  common  with  those  of  the  modern  tech¬ 
niques,  a  point  to  which  we  shall  return  later. 

There  seem  to  be  two  principal  reasons  why  acoustics  has  suffered  in  this  respect,  despite 
having  been,  a  hundred  years  ago,  the  field  in  which  MAE  almost  started.  These  are  firstly, 
that  classical  acoustics  appeared  by  around  1910  to  be  a  fairly  well  worked  out  subject,  largely 
because  of  the  great  work  done  in  the  second  half  of  the  1 9th  Century  by  Rayleigh  and  others. 
There  were,  of  course,  difficult  diffraction  problems  (such  as  the  quarter-plane  diffraction 
problem)  for  which  the  mathematical  techniques  of  1910  were  quite  inadequate  and  for  which 
even  today  they  are  hardly  adequate.  But  the  physical  principles  of  classical  acoustics  seemed 
well  enough  understood,  and  as  a  result  research  in  acoustics  was  taken  by  physicists  into 
directions  ever  more  remote  from  the  studies  of  the  compressive  wave  behavior,  at  reasonable 
frequencies,  of  air  and  water,  which  had  been  the  previous  concern  of  acoustics.  Now  of 
course  classical  acoustics  was  not  worked  out  by  1910,  but  it  took  another  forty  years  until 
Lighthill's  work  on  aerodynamic  noise  in  the  early  1950’s  for  it  to  be  realized  that  the  most 
essential  aspect  of  ordinary  acoustics  (the  aspect  of  energy  transfer  from  nonpropagating 
modes  of  one  kind  or  another  into  propagating  acoustic  energy)  was  in  fact  missing. 


I*  *  Technical  References  are  listed  at  the  end  of  the  text. 


That  aspect  now  pervades  most  work  on  ordinary  acoustics,  either  explicitly  or  implicitly,  and 
has  been  responsible  for  much  of  the  great  development  which  has  taken  place  in  acoustics  in 
the  last  25  years.  Nonetheless,  much  of  acoustics  is  still  in  the  hands  of  physicists  working  in 
very  different  areas,  and  physicists  have  been  conspicuously  slow  to  make  any  use  of  the 
perturbation  techniques  which  have  been  devised  by  engineers  and  applied  mathematicians 
working  in  fluid  dynamics. 

A  second  group  of  workers  has  also  been  slow  to  use  these  techniques.  Those  are  the 
workers  in  wave  diffraction  theory,  a  field  in  which  mathematical  rigor  is  often  necessary  and 
has  in  any  case  become  customary.  Against  the  usual  background  in  diffraction  theory  (of 
proofs  of  the  convergence  of  iterative  solutions  to  integral  equations,  of  the  solution  of 
problems  with  the  aid  of  function-analytic  arguments  in  one  or  more  complex  variable  planes, 
etc.),  techniques  like  MAE  appear  very  suspect  indeed.  There  is  not  much  that  can  be  proven 
at  intermediate  stages  of  the  calculations  and  the  validity  of  the  final  results  can  usually  only 
be  proven  in  rather  simple  cases.  If  one  is  not  worried  unduly  by  lack  of  mathematical  rigor, 
however,  the  gains  from  using  modern  perturbation  methods  can  be  immense.  I  quote  three 
examples,  in  the  first  of  which  at  three  successive  annual  Theoretical  Mechanics  Colloquia  in 
the  U.K.  were  presented  (i)  a  rigorous  derivation  of  two  terms  for  the  scattering  cross-section 
in  a  high-frequency  water  wave  problem,  running  to  50  pages  of  print:  (ii)  an  equally  rigorous 
derivation  of  the  same  two  terms  by  a  more  efficient  method  requiring  24  pages  of  print; 

(iii)  use  of  MAE  to  find,  in  2  pages  of  print,  four  terms  of  the  expansion  for  the  cross-section 
and  to  detect  an  error  of  sign  in  the  second  term  of  (i)  and  (ii).  As  the  second  example. 

Dr.  Frank  Leppington  and  1  have  used  MAE,  an  approximate  method,  to  derive  exact  closed 
form  solutions  for  certain  quantities  in  a  thick  plate  diffraction  problem  [  7 J  which  previously 
had  only  been  given  approximately  using  a  so-called  exact  modification  of  the  Wiener-Hopf 
technique.  As  the  third,  asymptotic  expansions  for  certain  modified  Mathieu  functions  were 
needed  18)  in  connection  with  a  study  of  the  behavior  of  a  fluid  jet  of  elliptic  cross-section  - 
modelling  the  kind  of  exhaust  jet  which  propels  the  CONCORDE  SST.  I  used  MAE  to  find 
these  expansions,  which  turned  out  to  be  new,  and  which  have  since  been  established  rigor¬ 
ously  by  one  of  my  colleagues,  Dr.  W.  Barrett,  of  Leeds  University,  England. 

I  hope  in  these  notes  (a)  to  convince  the  reader  that  singular  perturbation  techniques, 
and  in  particular  MAE,  form  a  particularly  valuable  tool  in  the  context  of  acoustic  problems 
and  (b)  to  show  through  simple  examples  in  acoustics  how  MAE  can  be  applied  in  a  reasonably 
systematic  and  safe  way.  Dr.  Martin  Lesser  and  I  have  attempted  this  elsewhere  [9]  in  an 
article  specifically  on  MAE  in  acoustics.  In  retrospect,  however,  that  article  seems  to  involve 
rather  too  much  formalism  in  its  attempt  to  make  the  method  routine,  and  seems  also  to  deal 


with  too  wide  a  range  of  problems  for  the  kind  of  course  needed  here.  Accordingly,  these 
notes  will  be  restricted  to  deal  with  a  few  examples  in  detail,  after  which  reference  to  the 
paper  (9)  may  be  valuable.  Our  main  aim  will  be  to  establish  firmly  the  basic  ideas  of  the 
MAE  technique  in  the  simplest  relevant  context. 


2.  HEURISTIC  APPROACH  TO  MATCHING  AND  ITS  PITFALLS 


A  heuristic  approach  to  matching  has  been  widely  used  in  the  past.  Our  purpose  here  is 
to  show  how  simple  acoustic  problems  demonstrate  that  anything  short  of  a  proper  mathemat¬ 
ical  rule  for  matching  is  likely  to  lead  to  erroneous  results,  even  it  the  physical  basis  for 
heuristic  matching  is  clear.  The  correct  solutions  will  be  found  later  with  the  aid  of  an 
Asymptotic  Matching  Principle. 

We  start  by  looking  at  the  solution  given  by  Landau  and  Lifschitz  [10)  for  the  sound 
field  radiated  into  still  fluid  by  the  forced  oscillation  of  a  body  about  a  fixed  mean  position. 
If  the  velocity  of  the  (rigid)  body  is  V  exp(-io)t),  we  want  a  solution  of 


with 

and 


(72  +  kg)0  -  0 

0  ~  r'1  exp(ikQr)  fid,  0)  as  t  -*<*> 
n  •  V0  =  n  •  V  on  the  body 


kQ  =  co/c0  being  the  acoustic  wavenumber  at  frequency  w.  If  the  body,  of  typical  dimension 
L  say,  is  compact- of  small  extent  compared  with  the  wavelength  27rk^'  -  then  we  argue  that, 
with  axes  fixed  on  an  origin  somewhere  in  the  body,  the  motion  is  incompressible  within  a 
wavelength  of  the  body.  Thus  for  r  «  k^1  we  have 


and 


720  =  0 
n  •  V0  =  n  •  V , 


(2.2) 


the  radiation  condition  being  inapplicable  for  this  range  of  r.  It  is  argued,  however,  that  0 
cannot  grow  when  r  55>  L  and  therefore  that  no  eigensolution,  which  would  necessarily  grow 
algebraically  for  r »  L,  can  be  present.  (The  term  eigensolution  is  used  for  any  solution  of 
V20  =  0  with  n  •  70  =  0  on  the  body.  For  a  sphere  of  radius  L,  for  example,  the  general 
axisymmetric  eigensolution  is 


»-£a„v«»*)W  nL--—  \ 

r^i  n  n  *  (n  +  l)  r"  +  1  ^ 
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where  Pn  denotes  the  nth  Legendre  polynomial.)  Then  for  r  »  L  this  inner  potential  4>  has 
an  harmonic  multipole  representation,  in  which  the  monopole  must  be  absent 


1  9Xj  \r/  yax^XjW 


(2.3) 


the  leading  term  of  which  is  a  potential  dipole, 

0~a‘£;(r)' 


(2.4) 


For  example,  in  the  case  of  a  sphere  of  radius  L  the  boundary  condition  on  <p  is 

30 

—  =  V  cos  6  on  r  =  L 
dr 


if  the  velocity  V  is  in  direction  0  =  0,  and  the  potential  is 


VL3 

<t>  = - cos  6 

2r2 


which  is  precisely  (2.3)  if  a  =  (1/2  VL3,  0,  0)  and  ay  = _ =  0. 

Now  we  argue  in  a  complementary  fashion  that  if  r  »  L  then  the  body  appears,  to 
leading  order,  as  a  singularity  at  the  origin,  so  that  for  r  »  L  we  seek  a  solution  of  the  full 
Helmholtz  equation,  satisfying  the  radiation  condition,  with  singular  behavior  at  r  =  0.  The 
most  general  such  solution  is 


<t>  =  b 


r  1  dXj 


+  b 


32 


y  dXjdxj 


(2.5) 


b  ,  .  ± 
r  ‘  3x: 


+  ... 


when  r  <5C  k^1 .  Thus  in  an  overlap  domain 


L  «  r  «  k’1 


appropriate  terms  from 


5 


should  be  identical.  The  obvious  choice  is 


b  =  0,  a;  =  bj ,  (2.6) 

and  thus,  since  a  is  known  in  terms  of  V  for  a  body  of  given  shape,  we  arrive  at  the  leading 
order  wave  field  in  the  form  of  an  acoustic  dipole 


~  ik 


o 


(a  •  x) 

— 7=-  exp(ikQr) 


(2.7) 


This  result  is  in  fact  correct  for  a  compact  body  of  any  shape.  The  argument  is  presented 
with  such  authority  by  Landau  and  Lifschitz  [10]  that  one  is  only  inclined  to  examine  its 
weak  points  when  one  uses  the  same  argument  only  to  arrive  at  a  plainly  wrong  result. 
Although  this  simple  kind  of  argument  was  used  a  lot  by  Lord  Rayleigh  [11]  and  by 
Lamb  (121  (to  examine  sound  transmission  through  rows  of  parallel  slits  in  a  screen,  for 
example),  it  is  not  difficult  to  find  simple  problems  in  which  a  strictly  comparable  result  goes 
badly  wrong,  and  we  now  give  two  such  examples. 

Take,  for  instance,  the  case  of  plane  wave  scattering  by  the  said  compact  body.  If  the 
incident  potential  is 

01  =  exp  ikQx 

then  the  scattered  field  <t>  satisfies 

(V2  +  k2)  <t>  *  0 

0  ~  r'1  exp(ikQr)  f{0)  as  r-»°°  (2.8) 

n  •  70  =  -  n  •  V01  on  the  body 

which  is  almost  the  problem  (2.1)  again.  For  r  «  k^1  the  potential  is  harmonic  and  the 
incident  stream  is  effectively  a  uniform  stream,  01  ~  (ik0)x  and  therefore  in  the  case  of  a 
sphere  the  leading  order  inner  field  is  precisely  as  before  if  we  make  the  substitution  V  =  -ik0. 
Whether  or  not  the  scatterer  is  a  sphere,  the  inner  potential  still  has  the  form  (2.3).  dominated 
for  r  »  L  by  an  harmonic  dipole.  Again,  we  can  evidently  “match”  this  dipole  to  an  outer 
wave-field  acoustic  dipole 


6 


.  a  (jkA 

;  3Xj  \  r  / 


by  choosing  the  strength  bj  of  the  acoustic  dipole  to  agree  with  the  aj  of  the  harmonic 
dipole  in  (2.4). 

This  time,  however,  the  casual  approach  is  quite  inadequate.  In  fact  the  directivity 
function  for  the  leading  order  wave  field  scattered  from  a  compact  sphere  is 

((d)  =  2-3  cos  d 

i.e.,  the  sum  of  equally  important  dipole  and  monopole  terms.  The  correct  monopole  strength 
cannot  be  found  from  first-order  matching  alone,  for  the  dipole  is  more  singular  than  the 
monopole  as  r-+0,  and  so  dominates  at  leading  order  as  we  come  in  towards  the  scatterer. 
The  monopole  strength  is  in  fact  determined  by  compressibility  effects  in  the  small  scattering 
region,  expressed  analytically  by  the  next  term  in  the  expression  near  the  body.  We  shall 
comment  on  this  again  after  looking  at  this  scattering  problem  in  detail  with  MAE. 

The  Landau-Lifschitz  problem  of  sound  emission  by  an  oscillating  body  also  goes  wrong 
when  there  is  a  uniform  mean  flow  past  the  body.  Then  the  answer  one  gets  depends  on 
whether  pressures,  or  potentials  -  or  some  other  field  variables  -  are  matched.  For  again,  the 
outer  behavior  of  the  potential  perturbation  near  the  body  is  a  dipole 


<t>  ~  a, 


and  the  leading  term  for  the  pressure  (p  =  -p(d/dt  +  U  3/dx)  <t>)  is 


(2.10) 


It  is  argued  that  this  should  be  matched  to  the  inner  limit  of  a  solution  to  the  convected 
wave  equation,  for  when  r  L  the  steady  perturbation  due  to  the  body  is  small  and  only 
the  uniform  flow  at  speed  U  remains.  For  simplicity  suppose  that  the  free  stream  Mach 
number  M  =  U/c0  is  small,  so  that  we  can  neglect  M2  compared  with  unity.  Then  a  wave- 
field  0  which  agrees  with  (2.9)  in  L  «.  r «  k'J  is 


(2.11) 


The  associated  pressure  in  the  wave  field  is  then 
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and  for  r  «  k'1  that  has  the  inner  limit 

p  ~ -p  aj(t)  ~  (1  -  M  cos  0)  (2.12) 

where  x  =  r  cos  0  and  the  stream  is  in  the  positive  x-direction.  Equation  (2.1  2)  does  not 
match  (2.10),  but  differs  from  it  by  a  Doppler  factor  (1  -  M  cos0).  Similar  differences  arise 
if  one  matches  the  pressures  and  then  considers  whether  the  velocities  agree. 

The  message  of  this  section  is  then  simply  this:  that  even  to  leading  order  some  care  is 
needed  in  ordering  terms  correctly,  and  a  mathematical  rule  for  matching  must  be  devised 
and  observed  without  the  use  of  physical  reasoning  to  shortcut  any  steps.  1  have  heard  it 
said  that  singular  perturbation  techniques  were  not  in  fact  invented  where  and  when  popularly 
supposed  (e.g.,  in  Cal  Tech  in  the  40’s  and  50’s)  but  by  such  men  as  Rayleigh  nearly  a 
hundred  years  ago.  1  tlunk  the  examples  given  here  show  that  that  is  not  a  tenable  view. 

The  great  19th  Century  physicists  discovered  a  physically  appealing  idea  with  something  in 
common  with  modern  methods;  but  they  did  not  discover  the  mathematical  formalism  for 
turning  their  ideas  into  a  reliable  method,  and  perhaps  did  not  appreciate  the  possibility  of 
errors  arising  in  the  way  1  have  tried  to  indicate.  It  is  regrettable  that  the  great  book  by 
Landau  and  Lifschitz  should  advocate  the  physical  rather  than  the  mathematical  approach, 
hut  it  is  interesting  to  note  that  they  attack  the  scattering  problem  by  a  completely  different 
method  (which  of  course  gives  the  right  answer). 

3.  FORMAL  APPROACH  TO  MATCHING 

It  is  a  fairly  general  criterion  that  a  perturbation  problem  (a  problem  containing  a  para¬ 
meter  e  where  we  are  interested  in  the  approximate  behavior  of  the  solution  as  e  tends  to 
some  value,  usually  zero  or  infinity)  is  singular  and  singular  in  a  way  calling  for  the  use  of 
MAE,  if  the  problem  involves  two  dynamically  significant  length  scales  whose  ratio  becomes 
either  large  or  small  as  e  approaches  its  limiting  value. 

Low-frequency  acoustic  generation  and  scattering  problems  exemplify  perfectly  this  kind 
of  singular  perturbation  problem  and  its  treatment  by  MAE.  (In  one  problem  1 7 1 ,  to  whi.  b 
we  shall  return  later,  an  acoustic  problem  served  to  highlight  an  extremely  subtle,  though 
probably  common,  failure  of  the  Asymptotic  Matching  Rule  as  widely  practiced  in  the 
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literature  and  to  show  how  a  modified  rule  could  be  devised  to  rectify  matters.)  The  pertur¬ 
bation  parameter  is  the  Helmholtz  number  e  *  k0L  which  tends  to  zero,  and  the  two 
obviously  significant  length  scales,  L  characterizing  the  body  geometry  and  k*1  characterizing 
the  wave  propagation,  become  asymptotically  disparate  as  e  -*  0. 

Consider  for  definiteness  j  scattering  problem,  in  which  r'  denotes  the  dimensional 
position  variable  and  ♦'(r’  k(>.  Lt  is  the  scattered  potential  for  an  incident  potential 
exptikox’>  Throughout  this  section  the  additional  dependence  on  angular  variables  (0,  0) 
will  be  suppressed  lo  assess  the  smallness  ot  various  terms  as  e  -*  0  we  need  dimensionless 
variables  I  he  potential  i-  dreads  dimensionless  and  either  L  or  k^1  can  be  used  to  normalize 
lengths  We  shall  use  the  following  notation 


r  k.  i  I-  •*  ‘o  'he  ■ u/er  position  coordinate 


l  “ 


♦<t  «  i  denotes  the  outer  potential 


R  -  -  -  -  denotes  the  inner  position  coordinate 
I  t 

♦'it  -  L.R.  k,.  L>  =  ♦tR.ti  denotes  the  inner  potential 

We  want  to  know  how  0(r.  el  behaves  as  e  -*  0  for  all  values  of  r  from  e  to  or 
equivalently  how  <I>(R  cl  behaves  as  e  -*  0  for  all  R  from  1  to  <»,  because  in  general  one 
needs  to  know  the  directivity  of  the  scattered  field  at  infinity  and  also  the  pressure  on  the 
scattering  body.  It  is  in  the  nature  of  a  singular  perturbation  problem  that  it  is  necessarily 
impossible  to  find  a  single  Poincare  asymptotic  expansion,  which  in  the  simplest  case  would 
be  a  power  series  in  e. 

0(r,  e)  ~  0Q(r)  +  e0,(r)  +  e2  02(r)  +  . . .  (3.1 ) 

which  holds  over  the  whole  range  of  values  of  r  of  interest.  This  is  because  when  such  a 
series  is  inserted  into  differential  equation  and  boundary  conditions  and  like  powers  of  e  are 
equated  one  has  to  assume  that  r  is  fixed  and  0(  1 ).  Usually  in  these  problems  r  can  be 
allowed  to  increase  indefinitely  without  affecting  the  process  of  equating  powers  of  e,  and  so 
boundary  conditions  at  infinity  -  the  radiation  condition,  for  example  -  can  be  imposed  on  the 
functions  0,0,,...  If  r  is  allowed  to  become  small,  however,  the  whole  basis  for  the 
expansion  (3. 1 )  is  undermined  as  the  terms  then  have  different  orders  of  magnitude  from 
those  postulated  in  the  expansion.  In  particular  this  means  that  the  terms  0O,  0,, . . .  cannot 
usually  be  made  to  satisfy  the  required  boundary  condition  on  the  body,  where  r  =  0(e). 


Consequently: 

(a)  We  cannot  find  the  pressure  on  the  scattering  body  from  the  outer  expansion 

0(r,  e)  ~  $Q(r)  +  e0, (r)  +  e2  02(r)  +  . . . 

which  is  unlikely  to  hold  when  r  is  as  small  as  e. 

(b)  We  cannot  at  the  moment  find  even  the  far-field  directivity  from  the  outer 
expansion,  because  that  expansion  contains  undetermined  constants,  and  even  functions, 
arising  from  the  fact  that  only  one  boundary  condition  (at  infinity)  can  be  applied  to  the 
outer  series. 

Similar  considerations  apply  to  the  inner  potential  which  also,  in  the  simplest  case, 
might  have  a  power  series  expansion 

4>(R,e)  ~  4>0(R)  +e<J>j(R)  +e2  <f>2 (R)  + . . .  (3.2) 

Formal  processes  (substitution  in  equations,  equating  powers  of  e,  etc.)  will  lead  to  a  partial 
determination  of  d>0,  4»j,  etc.,  only  for  values  of  R  of  order  unity.  One  expects  therefore  to 
be  able  to  apply  the  boundary  condition  on  the  body,  but  not  the  condition  at  infinity,  so 
that  we  can  certainly  not  find  the  directivity  function  from  the  inner  expansion  nor  can  we 
at  the  moment  even  find  the  pressure  on  the  scatterer  from  that  inner  expansion  because  of 
the  presence  of  undetermined  constants  and  functions  in  4>0,  <f>j,  . . . 

The  inner  and  outer  expansions  are  not  uniformly  valid.  The  outer  holds  in  the  wave- 
field  r  =  0(1)  and  usually  right  out  to  infinity.  It  may  also  hold  for  some  smaller  values  of  r, 
down  to  r  =  0(e1/2)  perhaps,  for  example,  but  it  will  not  usually  hold  uniformly  down  to 
values  as  small  as  r  =  0(e).  The  inner  holds  around  the  body,  R  =  0(1),  and  perhaps  also  for 
larger  values  of  R,  up  to  R  =  0(e'1/4)  say,  but  is  unlikely  to  hold  for  values  as  large  as  0(e*' ) 
which  are  in  the  wave  zone  where  r  is  the  natural  variable. 

Now  the  inner  and  outer  series  are  just  different  approximate  representations  for  the 
same  function  0'(r',  kQ,  L),  so  the  question  is  “Can  that  fact  be  utilized  to  pin  down  the 
indeterminacies  which  exist  in  the  inner  or  outer  expansions  separately,  and  so  to  find  approxi¬ 
mate  solutions  covering  the  whole  range  of  interest?” 

An  affirmative  answer  to  that  question  can  be  given  for  the  case  when  the  inner  and 
outer  expansions  overlap  strongly,  that  is,  when  both  are  simultaneously  valid  for  a  range  of 
intermediate  values  of  r,  r  =  0(e“)  say,  where 

0  <a,  <a  <a2  <  1  (3.3) 

The  number  a2  is  connected  with  the  smallest  value  of  r  for  which  the  outer  expansion  holds. 
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while  aj  is  connected  with  the  largest  value  of  r  for  which  the  inner  expansion  holds.  In  the 
example  just  quoted,  the  outer  expansion  held  down  to  r  =  0(eI/2>,  so  that  ctj  =  1/2,  the 
inner  up  to  R  =  0(e'l/4),  this  corresponding  to  r  =  0(e3^4)  and  giving  a,  =  3/4.  Since  aj  >or2 
there  is  no  overlap  in  this  case,  there  is  no  known  way  of  “matching”  the  inner  and  outer 
expansions,  and  indeed  no  general  method  may  exist  at  all.  If  the  inner  holds  up  to 
R  =  0(e"^2),  i.e.,  to  r  =  0(e^2),  then  Oj  =  ck2  =  1/2  and  there  is  only  marginal  overlap  along 
the  line  a  =  1/2  rather  than  in  a  domain  cXj  <  a  <  a2.  In  this  case  it  may  or  may  not  be 
possible  to  match  the  expansions.  No  general  rule  is  known  and  examples  can  be  given  of 
both  possibilities  [7,  1 3 j .  If,  on  the  other  hand,  the  outer  holds  down  to  r  =  0(e5/8)  say, 
the  inner  up  to  R  =  0(e'^2)  then  there  is  strong  overlap  in  the  domain  r  =  0(ea)  where 

1/2  <a  <  5/8 

Then  there  is  a  matching  rule,  which,  in  the  simpler  problems  at  any  rate,  is  sufficient  to 
uniquely  determine  the  various  unknowns  in  the  inner  and  outer  expansions. 

One  version  of  the  rule  runs  as  follows  (the  Intermediate  Matching  Principle).  Suppose 
one  knows  in  advance  that  an  overlap  domain  r  =  0(ea)  exists  for  a  given  by  (3.3).  (We 
should  note  that  in  general  the  values  oij  and  a2  will  depend  upon  the  orders  to  which  the 
inner  and  outer  expansions,  respectively,  are  carried;  a,  will  be  a  non-decreasing  function  of 
the  number  of  terms  retained  in  the  inner  expansion,  a2  will  be  a  non-increasing  function  of 
the  number  retained  in  the  outer  series,  so  that  the  width  of  the  overlap  domain  will,  if  any¬ 
thing,  reduce  as  the  expansions  are  taken  to  higher  order.)  Then  one  introduces  an  interme¬ 
diate  variable  r*  by  r  =  ear*.  R  =  r*/eua,  so  that  r  -*■  0  and  R  -*•  °°  as  e  -*•  0  with  r,  fixed. 
Inner  and  outer  series  are  then  expanded  as  e  -*■  0  for  fixed  values  of  r,  and  corresponding 
terms,  functions  of  r*  and  e,  are  made  to  agree.  If  one  does  not  know  the  range  of  permis¬ 
sible  values  of  a  in  advance  one  has  to  proceed  tentatively,  introducing  the  intermediate 
variable  and  trying  to  see  whether  there  is  a  range  of  values  for  a  which  will  allow  the  inter¬ 
mediate  expansions  of  the  inner  and  outer  series  to  be  matched  term  by  term. 

This  is,  at  present,  the  most  fundamental  of  the  available  matching  procedures,  though 
probably  not  the  most  common  nor  the  simplest  to  use  in  practice.  Difficulties  arise,  for 
instance,  when  the  inner  and  outer  expansions  are  re-expanded  in  terms  of  the  intermediate 
variable  r*,  it  not  being  clear  how  many  terms  in  the  expansion  are  to  be  retained  unless 
information  is  to  hand  on  the  domain  of  validity  of  the  inner  and  outer  series.  The  book  by 
Cole  {21  gives  an  extensive  illustration  of  the  use  of  the  intermediate  matching  idea,  although 
it  avoids  the  difficulties  inherent  in  this  approach  when  carried  to  high  order. 

The  most  common  procedure,  and  in  my  experience  by  far  the  most  straightforward  to 
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use  in  a  routine  fashion,  is  a  version  of  the  principle  put  forward  by  Van  Dyke  [  1  ]  as  the 
Asymptotic  Matching  Principle.  In  its  original  version  it  is  ambiguous  in  some  cases,  so  we 
propose  here,  as  in  (7),  to  use  the  following  notation  and  principles: 

We  write 

0(m)(r,  e)  =  0o(r)  +  e0,  (r)  + . . . .  +  em0m(r)  (3.4) 

for  the  outer  expansion  truncated  beyond  0(em).  If  the  outer  series  really  is  asymptotic,  then 

|0(r,e)-0<m>(r,e)|  =  o(em)  (3.5) 

for  appropriate  values  of  r.  Now  we  write  r  =  eR  in  hold  R  fixed  and  expand  through 
terms  0(en)  say.  This  gives  a  perfectly  definite  and  readily  calculated  set  of  terms  which  will 
be  denoted  by 

0(m*n)  (3.6) 

Next  we  take  the  inner  expansion  and  truncate  it  beyond  0(en),  giving 

«*»(") (R,  =  *0(R)  +  . . .  +  en4>n(R) .  (3.7) 

In  this  we  make  the  transformation  R  =  r/e,  keep  r  fixed  and  expand  as  e  0  through  terms 
0(em),  giving  another  set  of  terms  denoted  by 

(3.8) 

As  things  now  stand  is  a  function  of  R,  e,  4>^n’ a  function  of  r,  e,  but  we  shall 

understand  that  both  are  expressed  without  any  further  expansion  or  approximation  in  terms 
of,  say,  r,  e. 

Then  the  Asymptotic  Matching  Principle  states  that 

0(m,  n)  =  <j>(n,  m)  (3.9) 

It  is  easy  to  believe  the  truth  of  the  principle  on  the  basis  that  each  side  of  (3.9)  is  an 
asymptotic  representation  for  the  potential  in  the  overlap  domain,  and  asymptotic  representa¬ 
tion  with  respect  to  a  given  set  of  gauge  functions  like  integral  powers  of  e  is  unique,  so  that 
the  two  sides  of  (3.9)  must  be  identical.  That  reasoning  is  generally  quite  false,  though  the 
conclusion  as  expressed  in  (3.9)  is  correct.  Although  0^  and  do  have  asymptotic 
significance  with  respect  to  0',  the  functions  4m' n)  and  *<n>  in  general  have  no  asymptotic 
significance  whatever.  This  point  is  perhaps  made  clearer  if  we  consider  unequal  values  of  m 
and  n,  say  m  >  n.  Then  it  appears  that  genuine  information  has  been  thrown  away  in 


forming  <p^m' n)  from  while  information  of  no  significance  has  been  retained  in  forming 
<t>(n’  m>  from  The  principle  (3.9)  must  not,  however,  be  regarded  as  saying  anything 
about  the  common  asymptotics  of  the  two  expansions;  it  merely  defines  two  functions 
associated  with  the  expansions  in  a  clearly  prescribed  way,  and  asserts  their  equality. 

As  far  as  justification  of  the  MAE  procedures  goes,  one  can  say  the  following.  It  is 
necessary  first  to  prove  that  the  formal  substitution  of  an  assumed  series  like  (3.1,  3.2)  into 
equations  and  boundary  conditions  will  produce  an  asymptotic  series  for  any  values  of  r,  R. 
This  is  a  difficulty  common  to  all  perturbation  methods,  and  to  many  other  processes  (solution 
of  differential  equations  by  series  and  contour  integrals,  for  example)  as  well,  and  is  not  a 
difficulty  peculiar  to  MAE.  Under  this  heading  one  needs  also  to  establish  the  range  of  values 
of  r,  R  for  which  the  outer  and  inner  series  are  asymptotically  valid,  and  thus  to  find  the 
overlap  domain,  if  it  exists.  Proofs  of  these  points  can  be  given  in  some  cases,  but  the  proofs 
fall  far  short  of  covering  the  kinds  of  problems  of  immediate  research  interest.  Granted  the 
resolution  of  these  difficulties,  however,  the  two  overlapping  expansions  must  match  according 
to  (3.9).  This  has  been  proven  |6.  7.  13)  under  various  different  assumptions  about  the  over¬ 
lap,  and  with  greater  generality  than  we  have  so  far  required.  The  proof  given  in  [71  involves 
showing  that  there  exist  numbers  (m,,  n,  )  such  that  all  the  terms  inside  the  (m,  n)  “block” 

(the  one  of  interest)  plus  some  other  terms  outside  the  (m,  n)  block  but  inside  the  (mr  n, ) 
block  together  constitute  an  asymptotic  representation  for  <p'  with  an  error  smaller  than  the 
smallest  term  in  the  (m,  n)  block  and  for  some  intermediate  range  of  values  of  r.  This  is  done 
for  both  <t>  and  «l>  and  the  two  asymptotic  representations  are  equated  term  by  term.  The  rule 
(3.9)  is  just  a  statement  of  this  term  by  term  equality  for  a  particularly  efficiently  calculated 
block  of  terms. 

The  Asymptotic  Matching  Principle  has  been  proven  also  under  the  following  more 
general  conditions,  which  are  sufficient  to  cover  most  applications.  Suppose  that  the  inner 
and  outer  expansions  proceed  in  fractional  or  integral  powers  of  e  multiplied  by  integral 
powers  of  8ne, 

0(r,  e)  ~  0o(r)  +  e8ne  0,(r)  +  e02(r)  +  e28n2e  03 (r) 

+  e28ne  04(r)  e20s(r)  +  e38n3e  06(r)  +  . . .  (3.10) 


<l>(R,e)  ~  el/2<l>0(R)  +  e0, (R)  +  e3/28ne  02(R)  +  63/2<I>j(R) 
+  e2 8ne  4>4(R)  +  e205 (R)  +  e5/2 8n2e  06(R)  +  ... 


(3.11) 


being  a  case  that  arises  in  diffraction  by  a  thick  rigid  plate  [7)  and  is  interesting  in  several 
respects.  Then  a  function  like  $1 m )  or  4>^n)  will  be  constructed  as  before,  by  terminating  the 
expansion  beyond  0(em)  or  0(en)  as  the  case  may  be,  where  now  m  and  n  can  be  any 
rational  numbers  which  may  or  may  not  appear  in  either  of  the  series.  In  truncating  the 
series,  all  logarithmic  terms  are  to  be  grouped  together  according  to  the  order  of  their 
algebraic  multiplier,  and  are  never  to  be  separated  out.  Thus,  with  the  above  forms  for 
example,  we  have 

0tn  =  0O  +efine  0,  +e02 

0<1/2)  =  0O  . 

(*)(!)  =  fl/2  (K  + 

c  I 

4,13/2)  =  e>/2  4>^  +  e4>i  +  e-V2  gne  <j^  +  e3/2 

Having  taken  any  two  rational  numbers  m  and  n,  and  formed  and  4>*n\  one  then  writes 
r  =  eR  in  and  expands  up  to  and  including  all  terms  which  are  not  smaller  than  en, 
writes  R  =  r/e  in  4>^  and  expands  keeping  all  terms  not  smaller  than  em .  giving  the  quantities 
0<m,  n)  an£j  <j>(n.  m)  /{gain,  algebraic  order  only  is  taken  into  account.  e^Gn^e  being  regarded 
as  0(e^ )  whatever  the  value  of  p.  Then  the  matching  principle 

0(m,  n)  =  ^(n.  m)  (3.12) 

has  been  established  also  |7,  13]  for  expansions  of  this  kind. 

Occasionally  it  is  necessary  to  generalize  this  rule  further,  as  for  example  in  the  case  of 
scattering  by  a  soft  body  |‘)|.  But  that  example  deserves  separate  study  for  other  reasons, 
and  for  most  problems  the  principle  (3.12)  is  adequate. 

4.  SOUND  GENERATION  BY  FORCED  OSCILLATIONS 

We  return  now  to  the  problem  of  sound  generation  by  the  forced  oscillation  of  a  rigid 
compact  body.  For  definiteness  we  take  a  sphere  of  radius  L  whose  velocity  in  the  positive 
x'-direction  ( 0  =  0)  is  V  exp(-icot),  so  that  we  need  an  axisymmetric  solution  of 


-►  OO 


(4.1) 


In  terms  of  the  outer  variable  r  =  kQr'  and  with  an  arbitrary  normalization  for  0',  say 
0'  =  VL0,  the  complete  problem  is 


(V2  +  1)  0  =  0 


0  ~  r'1  exp(ir)  g(0)  as  r  -►  <» 


30  1  a 

—  =  -  cos  6  on  r  =  e 
dr  e 


} 


(4.2) 


with  e  =  k0L.  Letting  e  -*  0  makes  the  body  boundary  condition  tell  us  no  more  than  that 
0  must  be  singular  in  some  fashion  as  r  -*■  0.  Whatever  the  form  taken  by  the  outer  expansion 
(and  that  cannot  usually  be  assumed  in  advance,  for  it  is  dictated  to  some  extent  by  informa¬ 
tion  from  the  inner  region)  all  terms  in  the  outer  expansion  must  therefore  be  radiating 
axisymmetric  solutions  of  the  Helmholtz  equation,  the  general  solution  of  which  is  a  multipole 
series. 


a 


32 

a‘i  3Xj3x- 


+  . . . . 


(4.3) 


A  little  thought  shows  that  the  most  general  axisymmetric  solution  must  be 


+  am  V?  (  — J  +  (higher  multipoles  than  octupoles) 


*'22  3x  J_  \  r 

Let  us  assume  therefore  that  the  outer  series  starts  off  with  a  term  of  algebraic  order,  e^0o 
say.  where  0O  has  the  above  form,  in  which  denotes  the  Laplace  operator  in  the  two 
transverse  coordinates, 

V2  =  —  +  — 

1  3y2  3z2 


(4.4) 


The  coefficients  a,  a,,  a,,,  etc.  will  then  be  given  a  superscript  (o)  to  indicate  that  they  refer 
to0o. 

We  can  also  simplify  (4.4)  somewhat,  using  the  fact  that  for  r  ^  0 
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and  therefore  the  terms  with  coefficients  a,,  and  a,22  can  be  assumed  to  be  already  included 
elsewhere. 

Turn  now  to  the  inner  problem,  obtained  from  (4.2)  by  writing  4>(R,  e)  for  0(r,  e)  and 
r  =  eR,  so  that 


( Vr  +  e2 )  4>  =  0 

4>  ~  R'1  exp(ieR)  h(0)  as  eR  -*  °° 


94> 

—  =  cos  0  on  R  =  1 
oK 


\  (4.5) 


the  radiation  condition  obviously  being  unenforceable  in  the  limit  e  -*■  0. 

The  body  boundary  condition  indicates  that  (with  the  chosen  normalization,  <p‘  =  VL4>) 
4*  is  of  order  unity  in  the  inner  region,  so  that  the  inner  series  should  take  the  form 


4>  ~  4>  +  (terms  which  vanish  with  e). 


(The  only  other  possibility  would  be  that  4*  starts  with  a  term  which  becomes  infinite  as 
t  -*•  0,  4*  ~  e'1  4>,  say.  Then  4*  ,  would  have  to  be  an  inner  eigensolution,  the  general 
axisymmetric  form  of  which  was  given  following  equation  (2.2).  If  one  persists  with  this  and 
attempts  to  match  with  the  outer  solution,  one  finds  of  course  that  <I»_ j  =  0.  The  reader 
should  verify  this  after  seeing  how  the  first  nontrivial  matching  is  carried  out.)  The  problem 
for  <t>0  is  then 


Vr  -  0 


cos  0  on  R  =  1 


J 


(4.6) 


whose  general  solution  is 


<j>  =  _ 

T  r\ 


2R2 


oo 

cos  0  +  A<,0)  P^(cos  0)  |  R** 


1 


(4.7) 


the  sum  of  a  particular  solution  which  decays  at  infinity,  plus  a  general  axisymmetric 
eigensolution  satisfying 

,  94* 

Vr  4>  =  0,  =  0  on  R  =  I 
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For  matching,  note  that  when  <t>^  =  is  expanded  with  r  =  eR,  the  leading  term 
will  involve  only  inverse  powers  of  R,  while  when  4>(°)  =  <t>0  is  expanded  with  R  =  r/e,  the 
leading  term  will  involve  the  positive  power  xv  unless  we  choose  =  0.  Thus  no  inner 
eigensolutions  are  needed  to  leading  order  (though  they  are  needed  at  higher  order  to  match 
the  ascending  powers  of  R  which  arise  from  the  expansion  of  the  phase  factors  exp(ieR)  in 


(4.4)). 


We  then  have 


*(0)  =  <t»o  = 


(*)- 


- cos  0 

2R2 


- cos  0 

2r* 


and  so 


<j)(o,  o)  _  ^(o.  1 )  =  0 
2 

4,(0, 2)  =  _  co%0  =  «j>(o,3)  -  4>(o,4) 

2r2 

On  the  other  hand,  for  the  leading  outer  solution  we  have 

0(X)  =  e\ 

and  0^(eR)  =  0(e^-^)  where  N  is  the  highest  multipole  order  present  in  (4.4)  (N  =  1 
monopole,  N  =  2  dipole,  n  =  3  quadrupole,  N  =  4  octupole,  etc.).  Suppose  then  that  X  were 
equal  to  0.  Then  #(°)(eR)  =  0(e*' )  at  the  very  least,  and  because  that  term  must  match  with 
4>(°.o)  =  q  tj,e  on|y  possibility  is  that  0O  =  0.  The  same  conclusion  is  reached  if  we  take 
X  =  1,  because  =  0.  We  therefore  try  taking  X  =  2-i.e.,  the  order  of  the  leading  outer 

term  is  that  of  the  first  nontrivial  term  in  the  outer  expansion  of  the  leading  inner  term  -  so 
that  we  have 

0(2)  =  e2<t>0 . 

Since  the  expansion  of  4>(°^(r/e)  contains  no  inverse  power  of  e  it  follows  that  only  the  values 
N  =  1  and  N  =  2  are  permitted  in  0O,  and  hence 


«.-Kr)  t)- 


TTTS5 


(eR)  = 

..(0)  ± 

Ji  ax 

(£) 

(2.  0)  _ 

a*o)  = 

.,(»>  2 

-  cos  6 

>  3x\r; 

1  R 

2 

(2,  1)  _ 

.  O) 

(2.  2)  = 

_e2  (0)  JL(fL 

'  3X  \  2 , 

)  +<P{2' 

l)  =  _e2 

From  the  matching 


we  thus  get 


0(2.o>  =  4,(0.  2> 

Jo)  =  1 

ai  -> 


(4.10) 


which  determines  the  dipole  strength,  but  not  the  monopole  a^°\  We  see.  however,  that  this 
could  be  found  from  0(2, "  in  (4.9)  by  matching  to  4>(1'2),  an(j  therefore,  prompted  by  the 
fact  that  the  inner  expansion  of  0<2)(eR)  contains  terms  0(1)  and  0(e).  we  try  to  continue 
the  inner  expansion  with 


«f>  ~  4>0  +  e<J>,  +  o(e) . 


(4.11) 


4>.  is  an  inner  eigensolution 


and  the  expansion  of 

(7) +,*>  (7) 

must  start  with  terms  0<  I )  and  0(e)  in  order  that  4>(1,21  can  be  matched  to  </>(2-  ’*  in  (4.9). 
It  follows  that  only  the  term  with  v  =  1  can  be  present  in  (4.12)  and  that  then 


(4.12) 


jJjd.  2)  =  _  £_  CQS  Q  +  .  1  CQS  Q 

,  1  I 


(4.13) 


This  is  identical  with 
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if  and  only  if 


0(2.  1)  =  €a(o)  1 - L  cos  $ 

R  2R2 


a(D  =  a(°)  =  0  (4.14) 

Thus  matching  of  the  second  order  inner  solution  to  the  first  order  outer  solution  shows  that 
the  second  order  inner  solution  is  identically  zero  in  this  case,  and  that  the  previously  undeter¬ 
mined  leading  order  outer  monopole  strength  is  also  zero.  We  have  recovered  the  Landau- 
Lifschitz  result  (2.7),  and  the  “obvious”  choice  (2.6)  is  seen  to  be  correct -but  for  a  far 
deeper  reason. 

At  this  stage  we  will  break  off  and  make  a  similar  examination  of  the  plane  wave 
scattering  problem,  where  we  shall  see  that  the  second  order  inner  problem  is  not  trivial,  as  it 
was  above,  and  this  has  a  profound  influence  on  the  leading  order  outer  field,  making  the  so- 
called  “obvious”  choice  (2.6)  of  zero  monopole  strength  not  only  not  obvious,  but  actually 
incorrect. 


5.  PLANE  WAVE  SCATTERING 


For  the  scattering  by  a  fixed  rigid  sphere  of  radius  L  of  a  plane  incident  wave  of  poten¬ 
tial  exp(ikQx'),  the  problem  for  the  scattered  field  <t>  is 


(V'2  +kj)  0'  =  0 


0'  ~  r'_l  exp(ikQr')  f(0)  as  r'  -♦  «> 

(Ub* 

jp-  =  -ik0  cos  0  exp(ie  cos  0)  on  r'  =  L 


(5.1) 


The  potentials  are  already  dimensionless  and  hence  the  problem  in  outer  variables  takes  the 
form 


(V2  + 1)0  =  0 

0  ~  r*1  exp(ir)  g(0)  as  r  -*•  <» 

00 

=  -i  cos  0  exp(ie  cos  0)  on  r  =  e 
or 


(5.2) 


while  in  inner  variables  we  have 
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A 


(Vr  +€2)  4>  =  o 


<t>  ~  R'1  exp(ieR)  h(0)  as  eR  -*■  o° 


—  =  -  ie  cos  6  exp(ie  cos  6)  on  R  =  1 

oR 


As  in  £4,  all  outer  solutions,  to  any  order  in  e,  are  outer  eigensolutions.  with  the  general 
form  (4.4),  but  the  order  X  of  the  leading  outer  potential,  0  ~  e^0o,  is  not  yet  known.  It 
can  be  argued,  as  before,  that  the  order  of  the  leading  inner  potential  is  that  distinguished 
order  (a  term  used  by  Cole  (21)  for  which  the  solution  is  not  just  an  inner  eigensolution. 


(b  -w  £<f)  + 

o 


where 


oo 

•I*  -  +  —  COS  9  +  y*  A^Pytcos  0 )  |rp  +  (  1  (5> 

•’  ?r-  1  Vp+1/r«»  +  i/ 

and  again,  as  in  £4,  all  the  aJ,°*  =  0.  We  then  find  that  X  has  the  value  3  this  time,  with 
monopoles  and  dipoles  permitted  in  <t>Q,  so  that 

0  ~  e30o  +  . . . . 

*<>  =  ><0)  (7)  ♦  »i01  i  (7)  15  ■ 

It  might  be  thought  that,  since  X  is  3  here  rather  than  2,  quadruples  should  also  be 
permitted  in  0O  this  time.  That  cannot  be,  however,  for  then  e30o(eR)  would  contain  an 
0(1)  term  from  the  quadrupole,  which  could  not  be  matched  to  the  inner  solution  which  is 
0(e). 

The  matching  rule 


determines  the  dipole  strength  as 


0<  3.0=  $0.3) 


a<°>  =  -  i 


but  again  fails  to  determine  the  monopole  strength.  We  therefore  move  on  to  the  second 
inner  problem,  the  inner  expansion  continuing  with 
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<J>  ~  f4>Q  +e24>l  +  . . . 


The  difference  between  the  generation  and  scattering  problems  is  apparent  now,  in  that  4>,  is 
not  an  eigensolution  but  a  solution  of 


= 0 


=  cos2  0  on  R  =  1 


Writing  the  boundary  condition  as 


9*i  2  1 

—  >  jPj(COS«)  +  J 


we  easily  find  a  particular  solution  in  the  form 


P,  (cos0)  I 

4)  =  -  -  — - _  — 

1  9  r3  3R 


to  which  is  to  be  added  the  general  eigen  solution 


Application  of  the  matching  rule 


0<3,  2)  =  fp(2,  3) 


now  determines  the  monopole  strength  as 


„<°>  =  -  - 


and  shows  that  Aj/1)  =  0,  so  that  the  inner  eigensolutions  play  no  part  in  the  second  term 


The  leading  order  wave-field  is  then 


from  which  the  directivity  function,  as  r  -*■  then  follows  at  once  as  proportional  to  2-3 
cos  0,  as  quoted  in  f2.  Note  that  the  monopole  term  in  the  far-field  is  directly  comparable 
with  the  dipole  term,  yet  it  arises  from  a  weak  monopole  -e2/3R  in  the  near-field  as 


1 


compared  with  the  dipole  ie  cos0/2R2.  Such  subtleties  are  the  essence  of  Acoustics,  even  in 
the  low  frequency  limit,  and  Acoustics  is  therefore  an  excellent  vehicle  for  exposing  aspects 
of  techniques  which  are  easily  overlooked  in  apparently  more  complicated  but  basically  less 
delicate  fields. 

The  above  procedure  may  seem  rather  laborious  but  we  have  gone  into  it  in  some  detail 
because  a  routine  framework  of  this  kind  is  indispensible  in  more  difficult  problems.  More¬ 
over.  with  practice  the  way  in  which  the  expansions  are  likely  to  proceed  can  often  be  antici¬ 
pated  (though  as  Van  Dyke  [  1 )  emphasizes,  many  workers  have  obtained  incorrect  results  by 
failing  to  suspect  the  presence  of  certain  kinds  of  terms,  usually  logarithmic  terms,  whereas 
adherence  to  the  procedure  used  here  will  greatly  reduce  that  risk)  and  it  not  always  necessary 
to  include  the  most  general  eigensolution  if  one  keeps  in  mind  the  forms  of  the  terms  which 
need  to  be  matched  (though  again,  failure  to  find  the  most  general  eigensolution,  particularly 
in  nonlinear  problems,  has  led  to  many  incorrect  results  in  the  past,  and  general  eigensolutions 
should  be  included  if  there  is  the  slightest  suspicion  of  difficulty). 


6.  HIGHER  APPROXIMATIONS 


We  continue  here  with  the  scattering  problem  to  obtain  one  more  term  of  each  series. 
The  inner  series  so  far  is 

(  i  )  /  P2  <  cos  0 )  ,  \ 

cosl,r — —  t3R)+oif!'  1611 

If  <h(2,(r/e)  is  expanded  for  e  -*  0,  the  first  terms  are  0(t3)  and  0(e5 ).  so  we  attempt  to 
continue  the  outer  expansion  in  the  form 

<t>  ~  e30o  +  e5  +  o(e5 )  (6.2) 


using  the  matching  rule 


<J>(2.  5)  =  0(5.2) 


(6.3) 


Since  4>  =  0(e),  the  expansion  of  <t>(e R)  cannot  contain  any  terms  larger  than  0(e),  and  so  0, 
can  contain  monopole,  dipole,  quadrupole  and  octupole  elements,  with  the  general  form 


(6.4) 


Then 
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i 


.  «#**-» 


$<S)  (eR)  =  e30o  (eR) +  c50,  (eR) 


and  expanding  gives 


0<5)(eR) 


f 2  r2  2r2 

(1 +ieR-LjS- +  ..)  .  a  (1+ieR-^y- 


eR 


2e 


2  3X 


-} 


,  S  /  M>  -t-jgR  •«••••)  ^  3 

f  r  eR  ,2  ax 


(1  +  ieR 
R 


.,(1)  ,  .,(1)  , 

+  illL  At.  ( I  +  ieR  +  •  •)  am  33  (l+ieR  +  ...) 


e3  ax2  R  e4  ax3 

Truncating  this  expansion  beyond  0(e2)  gives 


R 


(6.5) 


(6.6) 


while  from  (6. 1 ) 


<J><2 


if3  1  P,  (cosff)  ,  I 

•  = -  COS  0  -  T.  e - Tf  — . 


2r2 


(6.7) 


Express  (6.6)  and  (6.7)  in  terms  of  the  same  variable  and  use  the  fact  that  P,  (cos  0)  = 
(1/2)  (3  cos2  0-1).  Then  the  matching  rule  (6.3)  requires 


—  cos0-^e2 
2R2  9 


1  ,  (3  cos2  0-1)  1  2  1  _  1  ,  1  ie 


R3 


I  i  2  1  -  1  2  1  j.  -  a 

'3e  R  =  ‘3e  R  +  ^C°S0 


9  cos  0  1  5  cos3  0~j 

R4  R4  J 

where  the  differentiations  in  (6.6)  have  been  performed  using  3R/3X  =  X/R  =  cos0.  Thus 
matching  gives 


+  e  a 


(n  (3  cos2  0-1)  (1) 


+  e  a 


RJ 


ill 


/ 


(6.8) 


23 


- —  — 


and  determines  the  most  singular  (octupole  and  quadrupole)  terms  in  ,  but  fails  to  deter¬ 
mine  the  dipole  and  monopole  terms. 

To  obtain  these  it  is  necessary  to  go  further  with  the  inner  expansion.  Consideration  of 
(6.5)  shows  that  the  dipole  term  in  (6.4),  with  coefficient  a  * 1  * ,  makes  a  contribution  to 
0<S)(eR)  of  order  0(e3),  while  the  monopole,  with  coefficient  a(,)  makes  a  contribution  of 
0(e4 ).  Therefore  a1,1*  will  be  determined  if  we  use  the  matching  rule 

0<s.  3)  _  $<s,  s>  (6.9) 


which  needs  the  third  term  of  the  inner  expansion,  while  the  monopole  can  only  be  determined 
from  the  fourth  term  of  the  inner  series  and  use  of  the  rule 

0<s.  4)  =  4,(4,  5)  (6 .10) 

Accordingly,  we  assume  an  inner  expansion 


which  gives  the  problem 


~  e  4>0  +  e2d>j  +  e34>2  +  e4«t3  + . . . 


VR  *2  =-*0=--2  COS* 


3<j> 

=  +  -j  cos3  6  on  R  =  1 


(6.11) 


(6.12) 


for  the  third  term.  This  is  typical  of  the  higher  order  inner  problems;  not  only  are  the 
boundary  conditions  inhomogeneous  (as  they  were  for  4>0  and  ),  but  the  Laplace  equation 
itself  now  has  a  forcing  function  related  to  the  lower  order  potentials.  Here  the  form  of  the 
forcing  suggests  looking  for  a  particular  integral  of  the  form  f(R)  cost),  and  then  one  quickly 
finds  that  f(R)  =  (i/4).  We  therefore  write 


and  then 


<*> 


2  =  —  COS  d  +  \j/ 


**o 


=  +1 

3R  2 


cos3  6 

P3(cos0)  +  ^  Pj(cosff) 


(6.13) 
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sag gateau 


(6.14) 


where  Pt  (cos  d)  =  cos  0,  P3  (cos  6)  =  (l/2)(5  cos3  6  -  3  cos  0)  and  the  transformation  to 
Legendre  functions  is  made  so  that  we  can  most  readily  use  the  facts  that  RnPn(cos0)  and 
R"(n  + 1 )  PR  (cos  d )  are  axisymmetric  solutions  of  Laplace’s  equation  (and  are  the  only  ones 
which  are  finite  on  the  axis,  where  cos  d  =  ±  I ).  We  now  try  to  eliminate  the  forcing  terms 
from  the  boundary  condition  (6.14)  by  looking  for  a  particular  solution 

P3(cos0)  Pj(cos0) 


which  gives  a  =  -  i/20,  0  =  -3i/20. 

To  this  we  have  to  add  the  general  inner  eigensolution,  as  given  in  (4.7)  with  coefficients 
A^,  so  that  the  general  solution  for  4>2  is 


:  :  P3 (cos  0) 

*2  =4°“  0-20—^— 


3^  Pt  (cos  6) 
20  R2 


Note  here  that,  whatever  the  values  taken  by  the  Ay2* ,  4*2  does  not  tend  to  zero  as  R  -*•  <» 
and  there  is  no  way  of  making  it  do  so.  At  still  higher  orders  the  inner  potentials  will  not 
merely  remain  finite  as  R  -»  »,  but  will  actually  tend  to  infinity. 

The  function  <J>(  3)  is  now  given  by 


4>( 


-0-S 


9  P,  (COS0)  3  .f3  ->■  P .(COS0) 

“s9-|‘5 -s — l;  +  Tcos9'2oe  ~~2 


i  Pj  (cos  6) 

20 C  14 


+  e 


in  which  the  leading  order  terms  must  be  at  most  0(e3)  in  order  for  this  to  be  matched  to 
the  outer  series  which  starts  with  e30o.  Therefore 

A<,2*=0  (6.17) 


and  again  it  seems  that  the  inner  eigensolutions  are  “too  singular’’  as  R  -*•  <*>  to  be  tolerated. 
We  shall  return  to  this  in  a  moment.  For  matching  according  to  (6.9)  we  have 
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*  • 


(6.18) 


i  „  2  ,P2(cos0>  f  iti  3j  P,  (cos  0) 

4A5.=e_c°s(i_5,!__ - _tTCOs9._e!__, 

-after  truncating  (6.16)  at  0(e5 )  and  then  returning  to  the  variable  R,  while  from  (6 .5)  we  have 
,<5  3)  1  e2  ie3  ic  P  cos  6  1  ,  .“1 

*  '-5  rTT  c  J 

<639, 

(use  having  been  made  of  the  results  (6  8)). 

The  equality  of  (6.18)  and  (6.19)  determines  the  dipole  coefficient  as 

a<,)  =  +  &,  (6.20) 


(6.20) 


and  all  other  terms  in  (6.19)  have  identical  counterparts  in  (6.18)  except  for  the  term  -ie3/3. 
This  outer  term  is  unmatched  because  of  a  failure  to  take  the  most  general  inner  eigensolution 
in  (6.15).  To  that  eigensolution  should  be  added  the  constant  term 


which  is  all  too  easily  discarded  as  irrelevant.  It  has  a  serious  effect  at  higher  order,  however, 
because  later  the  function  <1>2  will  act  as  a  forcing  function  for  the  inner  term  e54>4,  and  the 
constant  part  A^2'  of  4»2  will  induce  a  variable  part  of  d>4  and  hence,  through  the  matchings, 
of  the  wave  field  too.  It  we  add  A(Q2)  to  (6.15),  (6.18)  acquires  a  term 

e3A<2> 

which  can  be  matched  to  the  corresponding  term  in  (6.19)  provided 

A*,2*  --3.  (6.21) 

Thus  the  lowest  order  (least  rapidly  growing)  inner  eigensolution  is  called  into  play  in  the 
third  order  inner  solution. 

To  find  <f>3  in  (6. 1 1 )  we  have  the  problem 
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y 2  d>  =  _  * 
V  R  ^3 


2  Pj  (cos  9 )  j 

’  RJ  +3R' 


3  4*  2  , 

-r— -  =  -  —  cos4  9  on  R  =  1 
oK  o 


(6.22) 


This  is  a  complicated  but  straightforward  problem  which  is  tackled  as  before.  First  we  find  a 
particular  integral  to  annihilate  the  forcing  in  the  differential  equation.  Recall  that 

S  3R  (R2  3r)  +  ^2  ihTfl  dd  (s,n  6  m) 

and  that 

si?  55  {sin6id  Pn(cose))+n(n+l)Pn(cos0)  =  O, 
and  then  it  is  not  hard  to  construct  a  particular  integral 


so  that  we  write 


Then 


Ir  _l 

6*  27 


P2 (cos  6) 
R 


P2(cos0) 

R 


+  * 


VR  ^  _  o 


H  =  _ic°s<0-I- lp2(cos0)onR  =  ! 

Rewriting  cos4  0  as 

Js  P4  (cos  0)  +  ^  P2  (cos  0)  +  J 

we  find 

4  P4(cos0)  5  P2(cos0)  j 
*  =  525  ^5  +  567  ^3  +  5R 


(6.23) 


(6.24) 


plus  a  general  inner  eigensolution.  Since  the  outer  expansion  begins  with  0(e3),  e44>3(r/e) 
cannot  contain  any  terms  larger  than  e3  and  therefore  the  most  general  eigensolution  possible 
at  this  stage  is 


2 


f  i 


Hence 


A<3)  ♦  A<»  P^cos^jR  +  I^}. 


I  1  P«  (cos  6)  A  P.  (cos  0)  c  P,(cos0)  I 
4>,  =  -  R  -  TT  -=-= -  +  ~  - -  +  ^ 


2  6  27 


567 


5R 


525  Rs 

+  aJ,3)  +  A(,3)  P,(cos0){r  +  ^  —  1  (6.25) 

l  2  r2  J 

and  we  now  have  <I>(4),  that  is,  the  inner  expansion  through  0(e4),  determined  apart  from 
A^3)  and  A(,3). 

Next  we  calculate  <t>(4,5),  and  transform  it  back  to  the  variable  R,  with  the  result 

/  |  i  P,  (cosfl)  I 

4,(4. 5)  =  4,0, 5)  (eqn  6.!  8)  +  e4  1 1  R  _  ±  j— - +  -L 

(6.26) 


+  Aq31  +  A<3)  R  cos  0  | 


while 


0<S,4)  =  0(5.  3)(eqn  6>|9)+e4^_£l  ±(R2)+e4^.)  _fla<J>  J*L  (R).  (6.27) 

o  IZ  O A  K  Z  11 

Performing  the  differentiations,  using  (6.8)  for  a(, { 3  and  matching  (6.26)  to  (6.27)  gives 


l 

6 


1  (3  cos2  0  -  1 )  1 

54 


t  R  -  7T  x  — -  +  tt:  +  A'3)  +  A(3)  R  cos  0 


R 


5R 


siR'Rcos*+4%  «  (•  ss*l) 

6  6  R  18  \R  R  / 


so  that 


Aq3'  =  o 


l(3)  = 


.(l)  - 
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135 


(6.28) 


Thus  when  four  inner  terms  and  two  outer  terms  are  matched  these  terms  are  all  uniquely 
determined.  Note  that  in  the  fourth  inner  term  an  eigensolution  is  present  which  grows 
linearly  as  R  -*■  °o.  In  still  higher  approximations  the  further  eigensolutions  are  also  needed  in 
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the  inner  series  in  order  to  match  the  positive  powers  of  R  which  arise  from  the  expansion  of 
the  phase  terms  exp(ieR)  in  the  outer  series. 

There  is  no  particular  merit,  nor  any  great  difficulty,  in  carrying  the  inner  and  outer 
expansions  further,  and  so  one  can  determine  the  near  and  far  fields  to  any  order  in  e,  though 
the  algebraic  complexity  increases  rapidly  with  the  order.  The  point  we  want  to  demonstrate 
here  is  that  adherence  to  the  formal  procedures  provides  a  straightforward  and  routine  way  of 
carrying  the  approximations  through  a  number  of  terms,  and  that  the  matching  rule  can  be 
used  to  indicate  how  the  expansions  should  proceed  and  what  form  of  eigensolution  is 
sufficient  at  each  stage. 


7.  TWO  DIMENSONAL  PROBLEMS:  LOGARITHMIC  GAUGE  FUNCTIONS 


Logarithmic  gauge  functions  arise  in  many  problems,  and  the  aim  of  this  section  is  to 
show  how  they  arise  and  how  they  can  be  handled  by  MAE  in  the  context  of  two-dimensional 
acoustic  scattering  problems. 

We  start  by  looking  at  the  scattering  o!  a  plane  wave  01  =  expti  k0x')  by  the  rigid  fixed 
cylinder  r'  =  L,  so  that  the  problem  for  the  scattered  field  0'  is 


(V'2  +  k  i)<t>'  =  0 


00' 

^7  =  -  i  k0  cos  6  expti  k0L  cos  d)  on  r'  =  L 
0'  ~  (r')"1/2  expti  kQ  r')  f(0)  as  r’  -*  <» 


>  (7.1’ 


the  last  of  these  being  the  two-dimensional  radiation  condition.  In  terms  of  outer  and  inner 
coordinates,  we  have 


and 


(V2  +  1)0  =  0 


90 

9r 

0 


-  i  cos  0  exptie  cos  6)  on  r  =  e 
~  r',/2  exp(ir)  g(0)  as  r-+» 


>  (7.2) 


J 
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.  p-^iwn 


(V*  +e2)<t>  =  0 


—  *  -  if  cos  0  exp(ie  cos  0)  on  R  =  1 

oK 

<t>  ~  R'1/2  exp(ie  R)  h(0)  as  eR  -*■  °° 

the  scattered  field  being,  of  course,  symmetric  about  the  line  0=0. 
An  inner  eigensolution  A,,(R)  of  integral  order  v  is  defined  as 


Ay(R)  =  (R"  +  R-")  cos  v  0 


for  v  =  0,  1 ,  2 . and  satisfies 


7r  Aj,(R)  =  0 


-rr-  =  0  on  R  =  1  . 
dK 

The  general  outer  eigensolution  can  be  represented  in  various  forms.  The  commonest  repre¬ 
sentation  would  be  in  terms  of  a  sum  of  Hankel  functions  of  the  first  kind  and  of  all  integral 
orders. 


oo 

^  a„  H*,n  (r)  cos  v  0 


the  typical  term  of  which  varies  with  r  like  r'1/2  exp(ir)  for  large  r  and  so  satisfies  the 
radiation  condition.  Here  we  stick  to  our  previous  type  of  representation  in  terms  of  multi- 
pole  derivatives  of  the  fundamental  solution  H^V)  of  the  Helmholtz  equation  in  two 
dimensions.  Thus  we  write  the  general  outer  eigensolution  as 

+ ai  b,  H",(rl + a4  4^;  H°"(r’ + -  (7-5 

(x, .  x, )  being  the  position  coordinates  in  a  plane  perpendicular  to  the  axis  of  the  cylinder. 
As  before,  the  solution  symmetric  about  the  x,-axis,  0  =  0,  is  of  the  form 

aHo  >(r)  +  ai  ^  *>)♦.„  77  H^tr)  +  ...  (7.6 

dx 

An  advantage  of  using  the  form  (7.6)  is  that  we  need  only  to  know  the  behaviour  of  the 
zeroth  order  Hankel  function  as  r  -+  0,  which  is 
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H<,U(r)  ~  [l  -  ^  r2  +  0(r4)J 


+  —  (Snr  +  7F  -  fin  2) 

It  c 


r2+0(r4)  (7.7) 


where  >E  =  0.5772  ...  is  the  Euler  constant. 

As  in  the  previous  sections,  the  first  inner  problem  determines  how  the  expansions  start 
off.  Motivated  by  the  boundary  condition  in  (7.3)  we  assume 


and  then 


so  that 


—  e  <f>0  +  . . . 


*0  =  0 


=  -  i  cos  6  on  R  =  1 


®.  ■  nr  •  £  *!•"'  v*> 


and  unless  all  the  are  zero,  the  leading  term  of  the  expansion  of  4»(l)(r/e)  will  contain 
terms  which  are  either  constant  or  which  grow  with  r.  Such  terms  cannot  be  matched  to  the 
leading  order  solution,  because  that  is  an  eigensolution  like  (7.6)  whose  leading  order  terms 
all  contain  either  fin  r  or  inverse  powers  of  r.  Therefore  we  take 


4,0)  ig2  cos  § 


which  shows  that  the  leading  order  term  in  the  outer  expansion  must  be  0(e2),  so  that 


<t>  ~  e  0O  +  ••• 


(7.10) 
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where  0Q  is  given  by  (7.6)  with  superscript  ^  on  the  coefficients. 

Because  the  leading  inner  term  is  0(e),  e20o(eR)  cannot  be  larger  than  0(e),  and  so 


K  *  +  a<o)  j-  H^'(r) 


(7.11) 


is  the  most  general  possibility,  consisting  of  an  acoustic  monopole  and  an  acoustic  dipole 
respectively.  Then 

e20o(eR)  =  a^0)  e2  ~  ^Cn  R  +  fine  +  7E  -  Cn  2  - 

+  0(e4  Cne,  e4 )  +  e  a\0^  —  -^7—  +  0(e3  fine,  e3 )  (7. 1 2) 

1  7 r  K 


which  gives 


Applying  the  matching  rule 


*(2,n  =  ,  „  (o)  2i  cosj) 
6  f  Jl  7T  R 


0(2.1)  =  4,0.2) 


a(o)  =  - 
J'  2 


(7.13) 


(7.14) 


The  situation  is  precisely  as  in  the  spherical  problem.  First  order  matching  determines 
the  dipole  coefficient,  but  does  not  determine  the  monopole  part  of  the  first  order  outer 
field.  To  find  that  we  have  to  take  careful  account  of  small  phase  variations  in  the  scattering 
region  which  are  responsible  for  an  acoustically  efficient  source  mechanism.  Analytically,  we 
expand  (7.12)  to  next  order  to  see  how  the  inner  series  should  progress  and  find  that  the  next 
terms  are  0(e2  Cne)  and  0(e2),  just  as  indicated  in  (7.12).  We  therefore  assume 


4>  ~  e  4>0  +  e2  Cne  4>(  +  e24>2  +  . . . 


and  we  shall  assume  the  matching  rule 


*(2,  2)  —  ,2_(o) 


^  ^Cne  +  Cn  R  +  7E  -  Cn  2  -  ^  cos  $  =  <I>(2, 2) 


(7.15) 


(7.16) 


where  we  understand  that  e2  Cne  <!>,  and  e24>2  are  to  be  taken  as  a  single  term  0(e2). 

4>,  is  of  course  an  inner  eigensolution,  and  e2<t>,(r/e)  can  be  no  larger  than  0(e2),  since 
that  is  the  order  of  the  outer  potential.  Thus 
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(7.17) 


4>  =  A(1) 
no 


is  the  only  possibility.  For  4>2  we  have  the  problem 


K  *2  =  0 


d<l>. 


(7.18) 


-rr  =  cos'2  0  on  R  =  1 
oK 


This  problem  introduces  the  two-dimensional  source  potential  fin  R.  Writing  the  boundary 
condition  as 


94>2  j 

TiTI*1  +cos26) 


we  look  for  a  solution 


=  a  fin  R  +  0 


cos  20 


we  find  a  =  1 12.  0  -  -  1  /4,  so  that 


.  1  _  1  COS  2  0  -  2) 

*!-2t’'R-4-jr  +  Ao "• 


(7.19) 


the  form  of  the  eigensolution  being  determined  by  the  same  argument  as  for  (  7.17).  Thus 
we  have 

cos  2  0 


*(2.(R,  =  t„e  Arn  ♦  ,2 1  i  tnR  -  i  +  a'J1}  . 


*■«  (i)  =  t  e>  line  A»>  +  e’  {i  8n  ,  -  i  line  } 


tei  A(J> 


4  rJ 


(7.20) 


which  gives 


1 


+  e2  |  finr  +  e2  A^2) 


(7.21) 


Matching  (7.21)  to  (7.16)  gives 


iMu. 


...... 
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A(1>  =  - 
A°  2 


a(o)  _  *_ 
3  4i 


)  (7.22) 


Ao2)  =  I  (7e  -  £n  2  "  f) 


which  yields  the  monopole  strength  a(°)  and  non-zero  eigenfunctions  in  the  second  and  third 


inner  terms  4>(  and  <$2 . 


If  we  now  go  back  to  <l>(2)(r/e),  as  given  by  (7.20),  we  see  that  since  =  1/2  it  con¬ 

tains  terms  0(e2)  which  are  matched  to  the  outer  solution,  and  then  continues  with  an  0(e4) 
quadrupole  term  cos20/r2.  This  suggests  that  we  continue  the  outer  series  in  the  form 


<t>  ~  +  ... 


(7.23) 


The  choice  of  02  here,  rather  than  0, ,  is  deliberate.  02  is  again  an  outer  eigensolution,  and 
as  in  the  sphere  problem  we  can  allow  up  to  three  derivatives  of  H^’lr)  in  02,  but  not  more, 
in  order  that  e402(eR)  be  no  larger,  as  e  -*■  o,  than  the  leading  inner  term,  which  is  0(e). 
Thus 


02  =  a<2)  H^V)  +  a<2)  ~  H^r)  +  a<2)  il  H^r)  +  a<2{  ~  H^(r) ,  (7.24) 


3x 


and  without  taking  the  inner  solution  further,  the  only  matching  rule  we  can  use  is  that 

<I>(2- 4>(eqn.  7.20)  =  0(4>2). 


The  right  hand  side  of  this  equation  can  be  calculated  directly,  with  some  effort,  from  the 
assumed  form  (7.24)  and  the  expansion  (7.7).  One  finds  then  that 


„(2)  = 


Mil 


=  o, 


-lia<2>  =  .l 


(7.25) 


Ml 


and  the  monopole  and  dipole  coefficients  remain  undetermined. 

Expansion  of  0(4)  =  (e20o  +  e402)(eR)  beyond  0(e2)  produces  terms  0(e3  Gne)  and 
0(e3),  so  that  the  inner  expansion  must  take  the  form 


4>  ~  e  <t>0  +  e2  Gne  4>,  +  e24>2  +  e3  Gne  <t>3  +  e34>4  +  . . . 


(7.26) 


:  } 
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The  solutions  for  <I>3  and  <t>4  are 

4>j  =  A<3)  +  A<3)(R  +  R-1)  cos  6 


$4  =  A<,4)  +  A<4)(R  +  R'1 )  cos  6 


-  —  R  fin  R  cos  d  -77 

2  24 


cos  3  6 
R3 


Zi cos  2 

8  R 


and  the  matching  rule  to  be  applied  next  is 


(7.27) 


$(3.4)  =  0(4.3)  (7.28) 

Most  of  the  details  of  this  matching  are  unimportant.  One  aspect  serves,  however,  as  a  vital 
warning.  We  see  from  (7.27)  that  in  4>(3,4)  there  will  be  a  term 

-  ~  e4r  (En  r  -  Ene)  cos  0  (7.29) 


in  terms  of  outer  variables,  and  to  this  order  there  are  no  other  inner  terms  of  this  form.  In 
the  outer  expansion  (7.23)  there  are  also  no  terms  of  the  required  form;  certainly  there  are 
terms  which  contain  Ene  when  expressed  in  inner  variables,  but  when  these  terms  are 
reexpressed  in  outer  variables  the  lne  disappears  from  them.  Therefore  the  outer  expansion 
as  it  stands  contains  no  term  which  can  possibly  match  (7.29)  according  to  the  rule  (7.28)- 
and  (7.29)  has  a  definite  non-zero  coefficient,  because  the  term  in  (7.27)  which  gives  rise  to 
it,  namely  -  i/2  RfinR  cos0,  is  a  particular  integral  for  the  equation 


VR  ^4  M 


i  cos  8 

R 


The  fault  lies  with  the  assumed  form  (7.23)  for  the  outer  expansion.  Although  the 
expansion  of  4>(2*(r/e)  produced  only  algebraic  terms  0(e2,e4),  indicating  the  form  (7.23) 
the  presence  of  ine  terms  in  the  inner  expansion  at  one  stage  must  be  taken  as  a  warning  that 
they  will  probably  arise  at  the  next  outer  stage.  Thus  in  place  of  (7.23)  we  should  have 
anticipated 

0  ~  e20o  +  e4  Ene  0,  +  e402  +  . . .  (7.30) 


and  this  choice  enables  all  the  matchings  to  be  effected  without  difficulty  -  though  the 
algebraic  details  are  horrendous.  A  further  point  to  be  noted  is  that  because  0,  necessarily 
contains  Hankel  functions  which  have  logarithmic  singularities,  the  expansion  of  e4  Ene  0, 
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will  involve  a  terme4  En2e-and  therefore  the  inner  expansion  presumably  takes  the  form 
<t>  ~  e  4>0  +  e2  Ene  +  e2<P2  +  e 3  Ene  <f>3  +  e3d>4 

+  e4  En2e  <t>5  +  e4  Ene  <P6  +  e4<t>7  +  . . .  (7.3 1 ) 

It  is  necessary  in  fact  to  determine  all  the  eight  functions  in  (7.31)  in  order  to  uniquely 
determine  the  three  functions  in  the  outer  wave  field  (7.30). 

Once  again  then,  low  frequency  acoustic  problems  bring  out  subleties  in  the  matching 
procedure.  Although  we  have  not  proved  it  here,  it  can  be  shown  that  failure  to  regard 
e3Ene  d>3  and  e3<t>4  in  (7.261  as  a  single  “block”  leads  to  an  incorrect  determination  of  some 
of  the  coefficients  in  the  outer  terms  <p{  and  <p2  (see  [9]  and  [7]).  We  have  also  seen  how 
the  presence  of  logarithmic  terms  may  be  indicated  by  systematic  use  of  the  matching 
principle,  and  that  once  they  have  arisen  their  presence  must  be  suspected  in  all  subsequent 
terms  of  both  expansions. 


8.  PURELY  LOGARITHMIC  GAUGE  FUNCTIONS:  SCATTERING  BY  SOFT  BODIES 


The  problem  of  plane  wave  scattering  by  a  soft  cylinder- on  which  the  total  potential  is 
zero -provides  an  excellent  illustration  of  the  way  in  which  purely  logarithmic  gauge  functions 
may  arise,  leading  to  series  which  are  useless  for  practical  purposes  because  of  the  slowness  of 
their  convergence.  We  show  how  naive  application  of  the  matching  principle  fails,  but  then 
show  how  modification  of  the  gauge  functions  together  with  an  insistence  on  all  matchings 
enables  the  slowly  convergent  series  to  be  “renormalised”  into  a  rapidly  convergent  series. 

If  the  cylinder  is  circular,  of  radius  L,  then  the  inner  and  outer  problems  for  the 
scattered  field  $  corresponding  to  an  incident  potential  exp(ikQx')  are,  respectively. 


(V2  +  e2)  d>  =  0 
4>  =  -  exp(ie  cost?)  on  R  =  l 
(V2  +  I )  <t>  =  0 

<t>  ~  r'1^2  exp(ir)  1(0)  as  r  -*•  «> 
The  inner  eigensolutions  this  time  satisfy 


} 


(8.1) 


} 


(8.2) 


! 
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<t>  =  0  on  R  =  1  and  are 


*-o. 


A0  =  £n  R 


A„  =  (Rv  -  R^)  cosi>0  (v  =  1,2,...) 


while  the  outer  eigensolutions  are,  as  before,  of  the  general  form 


(8.3) 


>  +  »,  Jx  +  »„  ~  H>,"(r>  ♦  ... 

dx i 


(8.4) 


The  boundary  condition  in  (8.1)  indicates  that  <l>  =  0(1)  as  e  -*■  o,  so  that  if 


*  ~  *0  +  ■  •  • 


then 


OO 

4*0  =  - 1  40)  mr> 


v  =  o 


(8.5) 


and  we  cannot  decide  which  eigensolutions  may  be  permitted  without  looking  next  at  the 
outer  field.  Suppose  that  0  starts  off  with  an  algebraic  gauge  function  0  ~  ea  0O  with  a  >  0 
(on  physical  grounds).  Then  the  leading  term  of  ea  0o(eR)  must  be  0(  1 )  as  e  -*■  o  in  order  to 
match  4>0.  The  only  possibility  is  that  a  is  an  integer  and  that  the  first  term  permitted  in 
(8.4)  for  0O  contains  a  derivatives.  But  then  that  term  would  contain  a  factor  R-01,  whereas 
the  leading  term  in  the  expansion  of  4>0(r/€)  contains  only  r^  or  Bnr  where  0  >  0.  Matching 
to  leading  order  is  therefore  impossible  if  0  starts  off  with  an  algebraic  gauge  function.  A 
little  trial  and  error,  based  on  the  fact  that  a  monopole  a^H^Nr)  is  certain  to  be  present  in 
0O,  suggests  that  instead  we  should  try 

<8-6) 

Then  only  the  monopole  can  be  present  in  0O,  for  the  presence  of  a  dipole,  a,  3/dx  H^’lr) 
would  make  (l/8ne)  0o(eR)  infinite  like  (e  £ne)‘\  rather  than  0(1)  which  it  needs  to  be. 

Thus 

0O  -  a<°>  H(0u(r)  (8.7) 


and  with  an  obvious  extension  of  our  previous  notation. 
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(8.8) 


*(*“*)  (eR)  =  a<°> 


(eR) 


a(°)  2i , 


~  - - £ne  to  leading  order. 

fine  it 

Hoping  that  we  can  extend  the  matching  principle  to  cover  this  sort  of  case  we  write 


4> 


(frb’°)  =2i 


,(o) 


and  try  matching  this  to  the  terms  up  to  0(1 /fine)  in 

4>(°)  (l)  =  _  i  +  (fin  r  -  fine)  +  |  ^  j  )  jcosi>0. 

This  will  only  be  possible  if  all  the  A^,0^  =  0  and  then 

4>(°’  6ne)  =  -  1  , 


which  gives  us 


,(«)  =  +  ^ 


(8.9) 


This  determines  the  leading  order  solutions  as 

*~£^e7H°)(r) 


(8.10) 


To  improve  upon  these  approximations  we  have  to  expand  (l/£ne)  0o(eR)  up  beyond  the 
0(1)  term,  for  expansion  of  4>  =  -  1  is  not  helpful.  This  gives 


(eR)  =  -  1 


~{«nR  +  7E-8n2-f}  +  0(^i) 


and  suggests 


<l>  —  1  +  - —  +  . . . 

£ne  1 


(8.11) 


4>,  is  an  inner  eigensolution,  and  some  thought  shows  that  the  algebraic  type  of  eigensolution 
cannot  yet  enter  into  the  solutions.  (Note  that  our  previous  arguments  as  to  how  each 
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expansion  determines  the  form  of  the  other  have  to  some  extent  broken  down  with  purely 
logarithmic  gauge  functions;  for  example,  the  outer  expansion  of  4>0  =  -  1  is  -  1,  but  this 
does  not  mean  here  that  the  outer  series  start  with  an  0(1)  term).  We  try 


(8.12) 


(8.13) 


Clearly.  (8.13)  and  (8.14)  do  not  match,  though  the  choice  =  -  1  does  at  least 
match  the  term  in  (Cnr)/(Cn€).  We  can  either  be  satisfied  with  the  matching  that  does  work, 
and  which  determines  AjJ1)  uniquely,  and  merely  ignore  the  failure  of  other  terms  to  match, 
or  we  can  think  about  trying  to  modify  the  expansions  a  little  in  order  that  matching  of  all 
terms  can  take  place  satisfactorily. 

This  difficulty  often  arises  in  the  case  of  purely  logarithmic  gauge  functions.  Fraenkel 
[131  has  discussed  it  at  some  length,  showing  how  the  matching  failure  can  be  traced  to  at 
best  "marginal  overlap”  of  the  inner  and  outer  expansions  when  inverse  powers  of  2ne  are 
used  as  gauge  functions.  He  also  shows  how,  in  principle,  the  terms  which  cannot  be  expected 
to  match,  can  be  traced  at  each  stage.  It  should,  therefore,  be  possible  to  persist  with  our 
expansions  in  inverse  powers  of  2ne,  and  to  correctly  determine  all  terms,  regardless  of  the 
matching  failure,  by  careful  attention  to  the  points  made  by  Fraenkel.  But  the  result  would 
be  useless  for  practical  purposes,  for  the  coefficients  would  all  be  of  order  one  and  an  e  of 
I  O'10  would  still  only  make  consecutive  terms  smaller  by  a  factor  of  10  or  so. 

We  recommend  the  unrelenting  insistence  on  the  matching  rule,  combined  with  a  slight 
change  to  the  gauge  functions,  as  a  systematic  and  effective  way  of  redeeming  the  situation 
on  all  counts.  As  a  first  step,  we  try  taking  (£ne  +  K)'1  as  the  basic  gauge  function,  in  place 
of  (?ne)_l ,  K  being  a  constant.  Assume 
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<h - 6  +  . . . 

v  (8ne  +  K)  v° 


d>  ***  —  i  +  ■  d>  + 

*  (8ne  +  K) 


(8.15) 


and  as  before,  appropriate  choices  for  the  eigenfunctions  are 


0O  =a<°>  H^n(r), 


4>,  =  A^>  8n  R  . 


(8.16) 


icne-Hi) 


(eR)  =  (IneTKj  { 1  +  “ (6n  R  +  Cne  +  K  +  >e  "  6n  2  ~  K)  +  • } 


so  that 


K  ’  8ne  +  k)  =  +  [,  +  2i  + 

v  (8ne  +  K)  [_  n  K 


7e  -  8n  2) 


(8.17) 


jGsrnc  ■  5^)  ..  ,  . y.>  + 


0  (8ne  +  K) 


Equations  (8.17)  and  (8.18)  can  now  be  matched,  and  give 


a(°>  =  +  n 


A‘n=-1 


as  before,  while  the  remaining  terms  can  be  matched  if  we  choose 


K  =  7F  -  8n  2  -  -j  . 


Now  we  have  the  solutions 


(8.18) 


(8.19) 


(8.20) 


(8ne +  K) 


J _  .  I!  u(  i 

+  K)  2  H°  r  • 


<l> - I  - 


(8ne  +  K) ' 


(8.21) 
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and  to  take  the  expansions  further  we  expand  $(e R)  beyond  (fine  +  K)'1.  We  find 


$(eR)  —  1  - 


which  indicates  an  inner  expansion 


8n  R 


(fine  +  K) 


0  (tneTEi  '  e ’ ) 


$  -  -  1 


fin  R  e2  , 

-  + -  d>  +  + 

(fine  +  K)  (8ne  +  K)  2  3 


(8.22) 


This  shows  the  value  of  enforcing  the  matching  rules  in  terms  of  slightly  strained  gauge 
functions.  Instead  of  having  the  slow  series  in  inverse  powers  of  8ne,  we  have  effectively  now 
summed  that  series,  or  at  any  rate  a  subset  of  that  series,  into  the  term  (8ne  +  K),  and  the 
next  term,  0(e2),  is  now  very  much  smaller  than  1  /(Cne)2 ,  as  it  would  otherwise  have  been. 
Note,  of  course,  that  now  we  are  out  of  the  infinite  sequence  of  purely  logarithmic  functions 
we  must  expect  to  have  to  treat 


(TneTio  *1  and  e2*J 

together  as  effectively  a  single  0(e2)  term  in  the  matching. 

Failure  of  matching  at  higher  orders  in  this  and  similar  problems  indicates  that  a  further 
straining  of  the  gauge  functions  is  necessary;  the  overlap  has  become  marginal,  or  perhaps  has 
disappeared  altogether,  at  higher  order  in  terms  of  powers  of  e  and  inverse  powers  of 
(8ne  +  K).  The  remedy  may  be  to  try 


8ne  +  K 


8ne  +  K  +  K,e 


or  more  generally  8ne  +  K  -*•  8ne  +  K  +  K(  (e)  where  K(  (e)  -*•  o  as  e  -*•  o  and  is  to  be 
determined  by  the  enforcement  of  all  matchings. 


9.  COMPOSITE  EXPANSIONS 

Occasionally  it  is  useful  to  be  able  to  combine  the  inner  and  outer  expansions  into  a 
single  smooth  approximation,  valid  everywhere  from  the  scattering  body  to  infinity  in  the 
scattering  type  of  problem,  for  example.  (More  often,  however,  one  needs  just  to  know  the 
directivity  pattern  of  the  field  at  infinity,  or  the  pressure  on  the  scatterer,  and  these  can  be 
obtained  directly  from  the  outer  expansion  alone,  or  the  inner  expansion  alone,  respectively.) 
There  are  two  common  ways  of  combining  the  inner  and  outer  series  -  the  methods  of 
ADDITIVE  and  MULTIPLICATIVE  COMPOSITION.  In  the  additive  method  it  is  customary 


to  add  0^m^(r)  say  to  4>^(R)  and  take  off  the  part  0^m»n)  (or  4>(n,m\  they  are  identical) 
which  expansions  and  4>^n*  have  in  common  so  that  it  is  not  counted  twice.  This  gives 
a  composite 

0mn  =  ^(m)  +  *(n)  "  n)  (9. 1 ) 


which  is  often  claimed  to  be  as  good  as  the  inner  or  outer  series  in  their  respective  domains, 
and  so,  if  the  expansions  overlap,  to  provide  a  single  smooth  function  equivalent  in  accuracy 
to  the  two  separate  expansions.  In  the  method  of  multiplicative  composition  one  defines 


0(m)  .  ^(n) 

0(m’  n) 


(9.2) 


which  again  is  claimed  to  provide  a  different,  though  asymptotically  equivalent,  single  smooth 
function,  good  through  0(em)  in  the  outer  region  and  through  0(en)  in  the  inner. 

As  Van  Dyke  points  out  in  the  Notes  to  the  revised  edition  of  his  book,  these  two 
composites  are  particular  instances  of  a  very  general  class  of  composites.  Let  F  be  a  general 
functional,  F"1  its  inverse.  Then  a  general  composite  formed  from  the  outer  expansion 
through  0(em)  and  the  inner  expansion  through  0(en)  is  defined  by 

tf’mn  =  F'X  1  +  F(*(">)  -  F(0<m’  n>)  j  (9.3) 

The  functional  F(0)  =  0  gives  the  additive  composition,  the  functional  F(0)  =  8n0  gives  the 
multiplicative  rule.  Other  functionals  F  give  different  composites,  as  when  F(0)  =  exp  0  for 
example,  which  gives 

0mn  =  Gn  |  exp  0^m^  +  exp  -  exp  0(m’n)  |  (9.4) 

Now  in  general  none  of  these  composites  has  the  accuracy  claimed  for  it.  This  follows 
from  the  fact  that  in  general  neither  0(m>n)  nor  <|>(n’m)  has  any  asymptotic  significance  with 
respect  to  0  or  4>,  as  stated  already  on  pages  1 1  and  1 2  in  connection  with  the  Asymptotic 
Matching  Principle.  The  only  thing  that  can  be  said  about  the  additive  or  multiplicative 
composites  is  that  if,  for  example,  0*^  (r,  R  =  r/e,  e)  is  expanded  for  fixed  r  and  e  -*■  o,  and 
if  that  expansion  is  truncated  beyond  0(em)  the  result  will  be  precisely  0^m\  the  outer 
expansion  through  0(em).  Likewise  if  0*^  (r  =  eR,  R,  e)  is  expanded  through  0(en)  for 
fixed  R,  the  result  will  be  the  inner  expansion  through  0(en)-and  the  same  is  true  of 
the  multiplicative  composite.  The  same  is  not  true,  however,  for  any  functional  F  other  than 
the  additive  one  (F  =  1)  or  the  multiplicative  one  (F  =  8n)  because  of  the  presence  of  extra 
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terms  which  will  be  thrown  up  when  one  has  to  expand  not  simply  ^n\r/e),  for  example, 
but  F(<t>(n)(r/e)). 

Therefore  a  composite  is  at  best  only  as  good  as  or  4^  when  it  has  been  expanded 
and  truncated,  and  otherwise  may  not  have  the  accuracy  of  in  the  outer  region  or  of 
4>^  in  the  inner  region,  while  its  accuracy  or  otherwise  in  between -say  for  r  =  0(e,/2)-is  in 
general  quite  unknown.  This  last  point  is  emphatically  made  by  Schneider  (14)  who  gives  an 
example  in  which  the  multiplicative  composite  is  hopeless  in  the  intermediate  range  of  r 
because  of  a  zero  in  the  denominator,  Van  Dyke  [1,  Notes  to  revised  edition)  gives 

examples  in  which  different  functionals  F  can  be  used  to  produce  widely  differing  accuracies 
of  the  composites  formed  according  to  (9.3). 

My  feeling  at  the  moment  then  is  that  there  is  no  known  rule  which  invariably  produces 
a  composite  as  good  as  the  inner  or  outer  expansions  taken  separately,  let  alone  in  between, 
and  therefore  that  composites  are  to  be  avoided  wherever  possible.  In  simple  cases,  and  often 
to  leading  order  only,  it  may  be  possible  to  form  a  composite  by  careful  inspection  of  the 
behaviour  of  the  inner  and  outer  expansions  in  their  own  domains  and  in  the  overlap  domain. 
If  that  can  be  done,  the  additive  rule  should  be  used,  as  the  most  fundamental,  though  the 
multiplicative  rule  usually  gives  much  neater  results. 

Just  to  show  how  these  composites  are  formed,  according  to  the  usual,  though  generally 
incorrect  prescription,  we  go  back  to  the  field  scattered  by  a  compact  rigid  sphere,  for  which 


3  /  1  e'r  i  9  e,r  \  ( 3) 

*  'I'jT-iSTr* 


4>  ~  e  — —  cos  9  ~  4>(1) 
2R2 


(9.5) 


to  leading  order.  We  have 


1)  =  $>(1.  3)  =  e  — —  cos  0  =  <|>(  *3  f 
2R2 

so  that  additive  and  multiplicative  composition  both  produce  the  same  composite  which  is 
actually  0(3)  itself.  In  other  words 

g3 /  lei:  i±iL\ 

\  3  r  "  2  dx  t  f 

is  a  uniformly  valid  description  of  the  whole  field,  to  leading  order.  However,  the  accuracy  of 
this  leading  order  approximation  is  not  uniform,  as  is  only  to  be  expected  in  a  singular 
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r 


perturbation  problem.  In  the  outer  field  the  next  term  is  0(e5)  while  0(3)  is  0(e3),  whereas 
in  the  inner  field  0(3)  is  0(e)  and  the  next  term  is  0(e2). 

To  improve  accuracy  in  the  near  field  we  can  try  using 


with 


and 


Then  the  additive  (3,  2)  composite  is 


3  |  1  eir  i  3  eir  | 

7*2/2 

P2  (cos  6)  \ 

€  l  3  r  2  3x  r  ] 

7€  1 9 

R3  * 

and  for  once  this  is  neater  than  the  multiplicative  composite.  No  significance  is  to  be 
attached  to  the  last  term  here  in  the  wave  zone  r  =  0(1),  for  there  it  is  0(e5)  and  the  outer 
expansion  itself  is  only  good  to  0(e3)  so  far.  The  inclusion  of  the  last  term  when  r  =  0(1) 
may  make  a  distinct  difference  to  the  value  of  0  for  moderately  small  values  of  e,  and  may 
apparently  improve  upon  the  asymptotic  representation  0(3).  In  general,  however,  such  an 
improvement  is  coincidental,  and  the  032  composite  has  no  asymptotic  correctness  until  it  is 
expanded  for  fixed  r  and  terms  smaller  than  0(e3)  thrown  away,  which  of  course  takes  us 
back  to  the  outer  expansion  0(3)  alone. 

10.  CONCLUSIONS 

These  notes  have  attempted  to  show  MAE  at  work  on  some  simple,  but  hopefully 
representative,  problems  of  classical  linear  acoustics.  Despite  their  apparent  physical  and 
mathematical  simplicity,  these  low  frequency  sound  generation  and  scattering  problems 
illustrate  very  effectively  a  number  of  subleties  which  have  often  gone  unnoticed  in  more 
complicated,  but  perhaps  less  delicate,  problems  of  fluid  dynamics.  They  show  how  careful 
application  of  the  method  leads  in  a  reasonably  straightforward  way  to  solutions  in  the  near 
and  far  fields  whose  accuracy  is  limited  only  by  the  length  of  the  algebraic  manipulations  that 
have  to  be  carried  out  at  high  order. 


3  /  1  eir  i  3  elr)  (3) 

{*3  7-2  a;  vr* 

4>  ~e  —  cosfl  -e2{|  ?2  =  4 >(2) , 

2R2  *9  R3  3Rf 


0(3,  2)  =  <j,(2,  3)  =  _  f2  _L  +  e  — L_  cos  0 
3R  2R2 
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Modern  acoustics,  and  expeciaily  underwater  acoustics,  deals  with  many  problems 
involving  one  or  more  small  parameters.  For  example,  the  fluid  is  generally  only  weakly 
compressible,  so  that  often  the  Helmholtz  numbers  kQL  may  be  small;  a  vibrating  surface 
coupled  to  an  acoustic  fluid  may  be  subject  to  high  or  low  fluid  loading;  the  acoustic  fluid 
may  be  in  bulk  motion  at  a  low  Mach  number;  the  phase  speed  of  an  elastic  surface  wave 
may  be  only  slightly  less  than  the  sonic  speed,  or  the  frequency  may  be  very  much  less  than  the 
coincidence  frequency,  etc.,  etc.  All  such  problems  can  be  regarded  as  governed  by  two  distinct 
length  scales,  typically  an  acoustic  length  scale  and  the  other  scale  determined  by  hydro- 
dynamic  or  elastic  effects  regardless  of  compressibility.  As  such,  we  must  expect  them  to  be 
singular  perturbation  problems,  and  must  look  to  MAE  as  a  likely  method  of  attack.  A 
variety  of  problems  in  both  linear  and  non-linear  acoustics  is  discussed  from  this  viewpoint 
in  the  article  [91,  while  still  more  recent  and  unpublished  work  will  deal  with  MAE  in 
problems  involving  coupled  wave-bearing  media. 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

1  DTNSRDC  REPORTS.  A  FORMAL  SERIES  CONTAIN  INF ORN1ATION  OF  PERMANENT  r 
NICAL  VALUE  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  ' DE \  Tl F ICATlON  PEGARCN 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT 

2  DEPARTMENTAL  REPORTS  A  SEMIFORMAt  SERIES  CONTAIN  INFORMATION  OF  A  PF  • 
INARY.  TEMPORARY.  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OP  SlGNIf  ioA 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3  TECHNICAL  MEMORANDA  AN  INFORMAL  SERIES  CONTAIN  TECHNICAL  DOCUMEN TA 
OF  LIMITED  USE  AND  INTEREST  THEY  ARE  PRIMARILY  GORGING  PAPERS  INTENDED  F 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THE  I'  TYPE  AND 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT  ANY  DISTRIBUTION  OUTSID 
MUST  Bfc  APPROVED  RY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CA 
BASIS 
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